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Cold Dark Matter (CDM) has become the standard modern theory of cosmological structure formation. Its 
predictions appear to be in good agreement with data on large scales, and it naturally accounts for many properties 
of galaxies. But despite its many successes, there has been concern about CDM on small scales because of the 
possible contradiction between the linearly rising rotation curves observed in some dark-matter-dominated galaxies 
vs. the 1/r density cusps at the centers of simulated CDM halos. Other CDM issues on small scales include the 
very large number of small satellite halos in simulations, far more than the number of small galaxies observed 
locally, and problems concerning the angular momentum of the baryons in dark matter halos. The latest data 
and simulations have lessened, although not entirely resolved, these concerns. Meanwhile, the main alternatives 
to CDM that have been considered to solve these problems, self-interacting dark matter (SIDM) and warm dark 
matter (WDM), have been found to have serious drawbacks. 



1. Introduction 

The universe on the largest scales can be de- 
scribed by three numbers: 

• H = lOOhkms -1 Mpc~\ the Hubble pa- 
rameter (expansion rate of the universe) at 
the present epoch, 

• f2 m = p/pc, the density of matter p in units 
of critical density p c = SH^SnG)^ 1 = 
2.78 x 10 n /i 2 M o Mpc" 3 , and 

• Q\ = MSHq)^ 1 , the corresponding quan- 
tity for the cosmological constant. 

The current values of these and other key pa- 
rameters are summarized in the Table below (for 
additional references and discussion see [Q). It 
remains to be seen whether the "dark energy" 
represented by the cosmological constant A is re- 
ally constant, or is perhaps instead a consequence 
of the dynamics of some fundamental field as in 
"quintessence" theories [||. 

Cold Dark Matter (CDM) assumes that the 
dark matter is mostly cold — i.e., with negligi- 
ble thermal velocities in the early universe, either 
because the dark matter particles are weakly in- 
teracting massive particles (WIMPs) with mass 
~ 10 2 GeV, or alternatively because they are 



produced without a thermal distribution of ve- 
locities, as is the case with axions. CDM also 
assumes that the fluctuations in the dark mat- 
ter are adiabatic and have a nearly Zel'dovich 
spectrum. Considering that the CDM model of 
structure formation in the universe was proposed 
almost twenty years ago its successes are 

nothing short of amazing. As I will discuss, the 
ACDM variant of CDM with ft m = 1 - fi A ~ 0.3 
appears to be in good agreement with the avail- 
able data on large scales. Issues that have arisen 
on smaller scales, such as the centers of dark mat- 
ter halos and the numbers of small satellites, have 
prompted people to propose a wide variety of 
alternatives to CDM, such as warm dark mat- 
ter (WDM) H and self-interacting dark matter 
(SIDM) H. It remains to be seen whether such 
alternative theories with extra parameters actu- 
ally turn out to be in better agreement with data. 
As I will discuss below, it now appears that WDM 
and SIDM are both probably ruled out, while the 
small-scale predictions of CDM may be in better 
agreement with the latest data than appeared to 
be the case as recently as a year ago. 

In the next section I will briefly review the cur- 
rent observations and the successes of ACDM on 
large scales, and then I will discuss the possible 
problems on small scales. 
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2. Cosmological Parameters and Observa- 
tions on Large Scales 

The table below^] summarizes the current ob- 
servational information about the cosmological 
parameters, with estimated la errors. The quan- 
tities in brackets have been deduced using at least 
some of the ACDM assumptions. Is is appar- 
ent that there is impressive agreement between 
the values of the parameters determined by var- 
ious methods, including those based on ACDM. 
In particular, (A) several different approaches all 
suggest that Q, m w 0.3; (B) the location of the 
first acoustic peak in the CMB angular anisotropy 
power spectrum, now very well determined inde- 
pendently by the BOOMERANG § and MAX- 
IMA1 §] balloon data and by the DASI inter- 
ferometer at the South Pole 0, implies that 
Cl m + J7a ~ 1; and (C) the data on supernovae 
of Type la (SNIa) at redshifts z = 0.4 - 1.2 from 
two independent groups imply that £Ia — § H ~ j • 
Any two of these three results then imply that 
fl\ w 0.7. The la errors in these determinations 
are about 0.1. 

Questions have been raised about the reliabil- 
ity of the high-redshift SNIa results, especially 
the possibilities that the SNIa properties at high 
redshift might not be sufficiently similar to those 
nearby to use them as standard candles, and that 
there might be "grey" dust (which would make 
the SNIa dimmer but not change their colors). 
Although the available evidence disfavors these 
possibilities^] additional observations are needed 
on SNIa at high redshift, both to control sys- 
tematic effects and to see whether the dark en- 
ergy is just a cosmological constant or is per- 
haps instead changing with redshift as expected 
in "quintessence" models ||. But it is important 
to appreciate that, independently of (C) SNIa, 
(A) cluster and other evidence for fl m w 0.3 
jlj, together with (B) ~ 1° CMB evidence for 
Q m + Q\ ks 1, imply that J7a ~ 0.7. 

All methods for determining the Hubble pa- 

1 Updating the one in my talk at DM2000 0. 
2 For example, SNIa at z = 1.2 and ~ 1.7 apparently have 
the brightness expected in a ACDM cosmology but are 
brighter than would be expected with grey dust, and the 
infrared brightness of a nearer SNIa is also inconsistent 
with grey dust |ll|. 



Table 1 

Cosmological Parameters [results assuming 
ACDM in brackets] 



rameter now give compatible results, confirming 
our confidence that this crucial parameter has 
now been measured robustly to a la accuracy of 
about 10%. The final result H from the Hub- 
ble Key Project on the Extragalactic Distance 
Scale is 72 ± 8km s" 1 Mpc -1 , or h = 0.72 ± 0.08, 
where the stated error is dominated by one sys- 
tematic uncertainty, the distance to the Large 
Magellanic Cloud (used to calibrate the Cepheid 
period- luminosity relationship). The most accu- 
rate of the direct methods for measuring distances 
d to distant objects, giving the Hubble parameter 
directly as Ho — d/v where the velocity is deter- 
mined by the redshift, are (1) time delays between 
luminosity variations in different gravitationally 
lensed images of distant quasars, giving h « 0.65, 
and (2) the Sunyaev-Zel'dovich effect (Compton 
scattering of the CMB by the hot electrons in 
clusters of galaxies), giving h « 0.63 |lj,|l^]. For 
the rest of this article, I will take h = 0.7 when- 
ever I need to use an explicit value, and express 
results in terms of hjo = i/o/70kms _1 Mpc -1 . 

If tttot = 1 and structure formed from adiabatic 
initial conditions as assumed in ACDM, CMB 
data imply to = 14.0 ±0.5 Gyr |]|. For a ACDM 
universe with Q m = (0.2)0.3(0.4,0.5), the expan- 
sion age is t = (15.0)13.47(12.41, H.61)/if x Gyr. 



H = 100 h km s" 1 Mpc" 1 , h = 0.7± 0.08 
to = 13 ± 2 Gyr (from globular clusters) 

= [14 ± 0.5 Gyr, ACDM + CMB] 
fl b = (0.040 ± 0.002)/if 2 (from D/H) 

> [0.035fey 2 from Lya forest opacity] 
il m — 0.33 ± 0.035 (from cluster baryons etc.) 

= [0.34 ± 0.1 from Lya forest P(k)\ 

= [0.4 ± 0.2 from cluster evolution] 

w |n A - | ± | from SN la 
Cl tot = 1.04 ± 0.05 (from CMB peak location) 
£!a = 0.73 ± 0.08 (from previous two lines) 

< 0.73 (2a) from radio QSO lensing 
f2„ a 0.001 (from SuperKamiokande data) 

< [0.05 in ACDM-type models] 
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For f2 m w 0.3 and h « 0.7, there is excellent 
agreement with the latest estimates of the ages 
of the oldest stars in the Milky Way, both (A) 
from the globular cluster Main Sequence turnoff 
luminosities giving 12 — 13 ± 2 Gyr, (B) us- 
ing the thorium and uranium radioactive decay 
chronometers for halo stars jl?]], giving 14 ± 3 
Gyr and 12.5 ± 3 Gyr, respectively, and (C) from 
white dwarf cooling time, giving 12.7±0.7 Gyr as 
the age of the globular cluster M4 |jl8). It is re- 
markable that these four different clocks all agree! 

The lower limit on the hot dark matter (i.e. 
neutrino) contribution to the cosmological den- 
sity comes from the Super-Kamiokande atmo- 
spheric neutrino data pj|po| . The latest upper 
limit is from the 2dF redshift survey galaxy power 
spectrum 

3. Further Successes of ACDM 

The ACDM cosmology correctly predicts the 
abundances of clusters nearby and at z si within 
the current uncertainties in the values of the pa- 
rameters. It is even consistent with P(k) from 
the Lya forest |^| and from CMB anisotropies. 
Low-f2,„ CDM predicts that the amplitude of 
the power spectrum P(k) is rather large for 
fcs0.02/i/Mpc -1 , i.e. on size scales larger (k 
smaller) than the peak in P(k). The largest-scale 
surveys, 2dF and SDSS, should be able to mea- 
sure P(k) on these scales and test this crucial pre- 
diction soon; preliminary results are encouraging 
(23). 

The hierarchical structure formation which is 
inherent in CDM already explains why most stars 
are in big galaxies like the Milky Way [Q : smaller 
galaxies merge to form these larger ones, but the 
gas in still larger structures takes too long to cool 
to form still larger galaxies, so these larger struc- 
tures — the largest bound systems in the universe 
— become groups and clusters instead of galaxies. 

What about the more detailed predictions of 
ACDM, for example on the spatial distribution 
of galaxies. On large scales, there appears to 
be a pretty good match. In order to investi- 
gate such questions quantitatively on the smaller 
scales where the best data is available it is es- 
sential to do N-body simulations, since the mass 



fluctuations Sp/p are nonlinear on the few-Mpc 
scales that are relevant. My colleagues and I were 
initially concerned that ACDM would fail this 
test 1 24 , since the dark matter power spectrum 
Pdm(k) in ACDM, and its Fourier transform the 
correlation function (,i m (r), are seriously in dis- 
agreement with the galaxy data P g (k) and £ g (r). 
One way of describing this is to say that scale- 
dependent antibiasing is required for ACDM to 
agree with observations. That is, the bias pa- 
rameter b(r) = [£g(r) /£dm( r )] 1 ' 2 1 which is about 
unity on large scales, must decrease to less than 
1/2 on scales of a few Mpc This was 

the opposite of what was expected: galaxies were 
generally thought to be more correlated than the 
dark matter on small scales. However, when it 
became possible to do simulations of sufficiently 
high resolution to identify the dark matter halos 
that would host visible galaxies |26|,^7]], it turned 
out that their correlation function is essentially 
identical with that of observed galaxies! 

Jim Peebles, who largely initiated the study of 
galaxy correlations and first showed that £ g (r) ss 
(r/ro)- 1 - 8 with r » 5ft- x Mpc @, thought that 
this simple power law must be telling us some- 
thing fundamental about cosmology. However, 
it now appears that the power law £ 9 arises be- 
cause of a coincidence - an interplay between the 
non-power-law ^ m (r) and the decreasing survival 
probability of dark matter halos in dense regions 
because of their destruction and merging. But the 
essential lesson is that ACDM correctly predicts 
the observed £, g (r). 

The same theory also predicts the number den- 
sity of galaxies. Using the observed correlations 
between galaxy luminosity and internal velocity, 
known as the Tully-Fisher and Faber-Jackson re- 
lations for spiral and elliptical galaxies respec- 
tively, it is possible to convert observed galaxy lu- 
minosity functions into approximate galaxy veloc- 
ity functions, which describe the number of galax- 
ies per unit volume as a function of their inter- 
nal velocity. The velocity function of dark matter 
halos is robustly predicted by N-body simulations 
for CDM-type theories, but to connect it with the 
observed internal velocities of bright galaxies it is 
necessary to correct for the infall of the baryons 
in these galaxies J2£ - 31 , which must have hap- 
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pened to create their bright centers and disks. 
When we did this it appeared that ACDM with 
= 0.3 predicts perhaps too many dark halos 
compared with the number of observed galaxies 
with internal rotation velocities V w 200km s -1 
. While the latest results from the big sur- 
veys now underway appear to be in better agree- 
ment with these ACDM predictions [ ^5p^ ], this 
is an important issue that is being investigated in 
detail [37]]. Questions concerning the luminosity 
function still need to be resolved [ pq ]. 

The problem just mentioned of accounting for 
baryonic infall is just one example of the hy- 
drodynamical phenomena that must be taken 
into account in order to make realistic predic- 
tions of galaxy properties in cosmological theo- 
ries. Unfortunately, the crucial processes of espe- 
cially star formation and supernova feedback are 
not yet well enough understood to allow reliable 
calculations. Therefore, rather than trying to 
understand galaxy formation from full-scale hy- 
drodynamic simulations (for example [p9[ ), more 
progress has been made via the simpler approach 
of semi-analytic modelling of galaxy formation 
(initiated by White and Frenk |^b]^2]|, reviewed 
and extended by Rachel Somerville and me p3[ ). 
The computational efficiency of SAMs permits 
detailed exploration of the effects of the cosmo- 
logical parameters, as well as the parameters that 
control star formation and supernova feedback. 
We have shown that both flat and open 
CDM-type models with fl m — 0.3 — 0.5 predict 
galaxy luminosity functions and Tully-Fisher re- 
lations that are in good agreement with obser- 
vations. Including the effects of (proto-)galaxy 
interactions at high redshift in SAMs allows us 
to account for the observed properties of high- 
redshift galaxies, but only for Vl m w 0.3 — 0.5 pi| . 
Models with fl m = 1 and realistic power spectra 
produce far too few galaxies at high redshift, es- 
sentially because of the fluctuation growth rate 
argument mentioned above. 

In order to tell whether ACDM accounts in de- 
tail for galaxy properties, it is essential to model 
the dark halos accurately. The Navarro-Frenk- 
White (NFW) [|| density profile p NF w{r) oc 
r _1 (r + r s )~ 2 is a good representation of typical 
dark matter halos of galactic mass, except possi- 



bly in their very centers (§4). Comparing simula- 
tions of the same halo with numbers of particles 
ranging from ~ 10 3 to ~ 10 6 , my colleagues and I 
have also shown |4(J that r s , the radius where the 
log-slope is -2, can be determined accurately for 
halos with as few as ~ 10 3 particles. Based on a 
study of thousands of halos at many redshifts in 
an Adaptive Refinement Tree (ART) @ simula- 
tion of the ACDM cosmology, we jig found that 
the concentration c V i r = Rv'a/i"s has a log- normal 
distribution, with ler A(logc v j r ) = 0.14 at a given 
mass p9| , p0[ . This scatter in concentration re- 
sults in a scatter in maximum rotation veloci- 
ties of AV-max/Vmax = 0.12; thus the distribution 
of halo concentrations has as large an effect on 
galaxy rotation curves shapes as the well-known 
log-normal distribution of halo spin parameters 
A. Frank van den Bosch Q] showed, based on a 
semi-analytic model for galaxy formation includ- 
ing the NFW profile and supernova feedback, that 
the spread in A mainly results in movement along 
the Tully-Fisher line, while the spread in concen- 
tration results in dispersion perpendicular to the 
Tully-Fisher relation. Remarkably, he found that 
the dispersion in ACDM halo concentrations pro- 
duces a Tully-Fisher dispersion that is consistent 
with the observed onc.[] 

4. Halo Centers 

Already in the early 1990s, high resolution sim- 
ulations of individual galaxy halos in CDM were 
finding p(r) ~ r~ Q with a ~ 1. This behavior 
implies that the rotation velocity at the centers 
of galaxies should increase as r 1 / 2 , but the data, 
especially that on dark-matter-dominated dwarf 
galaxies, instead showed a linear increase with ra- 
dius, corresponding to roughly constant density 
in the centers of galaxies. This disagreement of 

3 Actually, this was the case with the dispersion in con- 
centration AHngf-.j ..) = 0.1 found for relaxed halos by 
Jing [kil, while we p8| found the larger dispersion men- 
tioned above. However Risa Wechsler, in her dissertation 
research with me feOJ , found that the dispersion in the 
concentration at fixed mass of the halos that have not had 
a major merger since redshift z = 2 (and could thus host 
a spiral galaxy) is consistent with that found by Jing. We 
also found that the median and dispersion of halo concen- 
tration as a function of mass and redshift are explained by 
the spread in halo mass accretion histories. 
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theory with data led to concern that CDM might 
be in serious trouble |^M^3[ . 

Subsequently, NFW |4q| found that halos in 
all variants of CDM are well fit by the pnfw{t) 
given above, while Moore's group proposed an al- 
ternative pm{t) oc r~ 3 / 2 (r + tm) -3 ^ 2 based on a 
small number of very-high-resolution simulations 
of individual halos p^-p^]. Klypin and collabo- 
rators (including me) initially claimed that typ- 
ical CDM halos have shallow inner profiles with 
a « 0.2 p9| , but we subsequently realized that 
the convergence tests that we had performed on 
these simulations were inadequate. We now have 
simulated a small number of galaxy-size halos 
with very high resolution (4(| , and find that they 
range between pnfw and pu- Actually, these 
two analytic density profiles are almost indis- 
tinguishable unless galaxies are probed at scales 
smaller than about 1 kpc. 

Meanwhile, the observational situation is im- 
proving. The rotation curves of dark matter dom- 
inated low surface brightness (LSB) galaxies were 
measured with radio telescopes during the 1990s, 
and the rotation velocity was typically found to 
rise linearly at their centers J57|-|5{J . But a group 
led by van den Bosch f60|| showed that in many 
cases the large beam size of the radio telescopes 
did not adequately resolve the inner parts of the 
rotation curves, and they concluded that after 
correcting for beam smearing the data are on the 
whole consistent with expectations from CDM. 
Similar conclusions were reached for dwarf galax- 
ies Swaters and collaborators showed that 
optical (Ha) rotation curves of some of the LSB 
galaxies rose significantly faster than the radio 
(HI) data on these same galaxies ]62|], and these 
rotation curves (except for F568-3) appear to be 
more consistent with NFW p3[ . 

Recently, a large set of high-resolution optical 
rotation curves has been analyzed for LSB galax- 
ies, including many new observations p4| . The 
first conclusion that I reach in looking at the den- 
sity profiles presented is that the NFW profile of- 
ten appears to be a good fit down to about 1 kpc. 
However, some of these galaxies appear to have 
shallower density profiles at smaller radii. Of the 
48 cases presented (representing 47 galaxies, since 
two different data sets are shown for F568-3) , in a 



quarter of the cases the data do not probe inside 
1 kpc, and in many of the remaining cases the res- 
olution is not really adequate for definite conclu- 
sions, or the interpretation is complicated by the 
fact that the galaxies are nearly edge-on. Of the 
dozen cases where the inner profile is adequately 
probed, about half appear to be roughly consis- 
tent with the cuspy NFW profile (with fit a aO.5), 
while half are shallower. This is not necessarily 
inconsistent with CDM, since observational biases 
such as seeing and slight misalignment of the slit 
lead to shallower profiles [ |65[ . Perhaps it is sig- 
nificant that the cases where the innermost data 
points have the smallest errors are cuspier. 

I think that this data set may be consistent 
with an inner density profile a ~ 1 but probably 
not steeper, so it is definitely inconsistent with 
the claims of the Moore group that a z 1.5. But 
recent work |36j has shown that Moore's simu- 
lations did not have adequate resolution to sup- 
port their claimed steep central cusp; the highest- 
resolution simulations appear to be consistent 
with NFW, or even shallower with a « 0.75. Fur- 
ther simulations and observations, including mea- 
surement of CO rotation curves |67[| , may help to 
clarify the nature of the dark matter. 

It is something of a scandal that, after all these 
years of simulating dark matter halos, we still 
do not have a quantative — or even a qualita- 
tive — theory explaining their radial density pro- 
files. In her dissertation research ]5C|] , Risa Wech- 
sler found that the central density profile and the 
value of r s are typically established during the 
early, rapidly merging phase of halo evolution, 
and that, during the usually slower mass accre- 
tion afterward, r s changes little (see also p8|). 
The mass added on the halo periphery increases 
R v ir-, and thus the concentration c V { r = R V i r /r s . 
Now we want to understand this analytically. 
Earlier attempts to model the result of sequences 
of mergers (e.g., |35^|7(J) led to density profiles 
that depend strongly on the power spectrum of 
initial fluctuations, in conflict with simulations 
(e.g. (f|]). Perhaps it will be possible to improve 
on the simple analytic model of mass loss due 
to tidal stripping during satellite inspiral that we 
presented in |j2|. Including the tidal puffing up 
of the inspiralling satellite before tidal stripping 
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can perhaps account for the origin of the cusp 
seen in dissipationless simulations, independent 
of the power spectrum 72 1 . They argue that the 
profile must be steeper than a = 1 as long as 
enough satellites make it into the halo inner re- 
gions, simply because for flatter profiles the tidal 
force causes dilation rather than stripping. The 
proper modeling of the puffing and stripping in 
the merger process of CDM halos may also pro- 
vide a theoretical framework for understanding 
the observed flat cores as a result of gas pro- 
cesses; reionization and feedback into the bary- 
onic component of small satellites would make 
their cores puff up before merging. This could 
cause them to be torn apart before they penetrate 
into the halo centers, and thus allow a < 1 cores 
|7l| . Other possible explanations for flatter cen- 
tral density profiles involving the baryonic com- 
ponent in galaxies has recently been proposed, in 
which the baryons form a bar that transfers angu- 
lar momentum into the inner parts of the halo[[74|, 
or alternatively binary black holes eject matter by 
a gravitational slingshot effect J75|. While these 
phenomena could be very important in massive 
galaxies, it is not clear that they are important 
in dark-matter-dominated dwarf and LSB galax- 
ies that have small or nonexistent bulge compo- 
nents. 

ft would be interesting to see whether CDM 
can give a consistent account of the distribution 
of matter near the centers of big galaxies, but 
this is not easy to test. One might think that big 
bright galaxies like the Milky Way could help to 
test the predicted CDM profile, but the centers of 
such galaxies are dominated by ordinary matter 
(stars) rather than dark matter 

5. Too Much Substructure? 

Another concern is that there are more dark 
halos in CDM simulations with circular veloc- 
ity V c £30 km s _1 than there are low-V c galaxies 
in the Local Group |78|]79f l . A natural solution 
to this problem was proposed by Bullock et al. 



4 Navarro and Steinmetz had claimed that the Milky Way 
is inconsistent with the NFW profile [79], but they have 
now shown that ACDM simulations with a proper fluctu- 
ation spectrum are actually consistent with the data [JFq]. 



||80|| , who pointed out that gas will not be able 
to cool in V c s30 km s _1 dark matter halos that 
collapse after the epoch of reionization, which oc- 
cured perhaps at redshift z re i on « 6 |pH - When 
this is taken into account, the predicted number 
of small satellite galaxies in the Local Group is 
in good agreement with observations (8^^|. ft 
is important to develop and test this idea fur- 
ther, and this is being done by James Bullock and 
by Rachel Somcrvillc and their collaborators; the 
results to date (e.g. |^3|,§4|) look rather promis- 
ing. Other groups (e.g. ]85|-|87j]) now agree that 
astrophysical effects will keep most of the sub- 
halos dark. As a result, theories such as warm 
dark matter (WDM), which solve the supposed 
problem of too many satellites by decreasing the 
amount of small scale power, may end up predict- 
ing too few satellites when reionization and other 
astrophysical effects are taken into account plfl . 

The fact that high-resolution CDM simulations 
of galaxy-mass halos are full of subhalos has also 
led to concerns that all this substructure could 
prevent the resulting astrophysical objects from 
looking like actual galaxies |79f| . In particular, 
it is known that interaction with massive satel- 
lites can thicken or damage the thin stellar disks 
that are characteristic of spiral galaxies, after the 
disks have formed by dissipative gas processes. 
However, detailed simulations |^|[k| have shown 
that simpler calculations had overestimated 
the extent to which small satellites could dam- 
age galactic disks. Only interaction with large 
satellites like the Large Magellanic Cloud could 
do serious damage. But the number of LMC-size 
and larger satellites is in good agreement with 
the number of predicted halos |78|, which sug- 
gests that preventing disk damage will not lead 
to a separate constraint on halo substructure. 

6. Angular Momentum Problems 

As part of James Bullock's dissertation re- 
search, we found that the distribution of specific 
angular momentum in dark matter halos has a 
universal profile fl92|| . But if the baryons have 
the same angular momentum distribution as the 
dark matter, this implies that there is too much 
baryonic material with low angular momentum 
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to form the observed rotationally supported ex- 
ponential disks (9^|^| . It has long been assumed 
(e.g. |p9 30|) that the baryons and dark matter 
in a halo start with a similar distribution, based 
on the idea that angular momentum arising from 
large-scale tidal torques will be similar across the 
entire halo. But as my colleagues and I argued re- 
cently, a key implication of our new picture of an- 
gular momentum growth by merging |94j ] is that 
the DM and baryons will get different angular 
momentum distributions. For example, the lower 
density gas will be stripped by pressure and tidal 
forces from infalling satellites, and in big mergers 
the gaseous disks will partly become tidal tails. 
Feedback is also likely to play an important role, 
and Mailer and Dekel ]7j| have shown using a 
simple model that this can account for data on 
the angular momentum distribution in low sur- 
face brightness galaxies 95 1. 

A related concern is that high-resolution hydro- 
dynamical simulations of galaxy formation lead to 
disks that are much too small, evidently because 
formation of baryonic substructure leads to too 
much transfer of angular momentum and energy 
from the baryons to the dark matter |)6). But 
if gas cooling is inhibited in the early universe, 
more realistic disks form |97| , more so in ACDM 
than in f2 m = 1 CDM |Q . Hydrodynamical sim- 
ulations also indicate that this disk angular mo- 
mentum problem may be resolved if small scale 
power is suppressed because the dark matter is 
warm rather than cold |9^], which I discuss next. 

7. Alternatives to ACDM? 

Because of the concerns just mentioned that 
CDM may predict higher densities and more sub- 
structure on small scales than is observed, many 
people have proposed alternatives to CDM. Two 
of these ideas that have been studied in the great- 
est detail are self-interacting dark matter (SIDM) 
§ and warm dark matter (WDM). 

Cold dark matter assumes that the dark matter 
particles have only weak interactions with each 
other and with other particles. SIDM assumes 
that the dark matter particles have strong elas- 
tic scattering cross sections, but negligible anni- 
hilation or dissipation. The hope was that SIDM 



might suppress the formation of the dense central 
regions of dark matter halos, although the large 
cross sections might also lead to high thermal con- 
ductivity which drains energy from halo centers 
and could lead to core collapse [ 100 1 , and which 
also causes evaporation of galaxy halos in clus- 
ters, resulting in violation of the observed "fun- 
damental plane" correlations [ 101 1 . But in any 
case, self-interaction cross sections large enough 
to have a significant effect on the centers of 
galaxy-mass halos will make the centers of galaxy 
clusters more spherical [102 103 and perhaps also 
less dense [104 ,[105 than gravitational lensing ob- 



servations | p_06| | indicate. 

Warm dark matter arises in particle physics 
theories in which the dark matter particles have 
relatively high thermal velocities, for example be- 
cause their mass is psI keV |107|, comparable to 
the temperature about a year after the Big Bang 
when the horizon first encompassed the amount of 
dark matter in a large galaxy. Such a velocity dis- 
tribution can suppress the formation of structure 
on small scales. Indeed, this leads to constraints 
on how low the WDM particle mass can be. From 
the requirement that there is enough small-scale 
power in the linear power spectrum to repro- 
duce the observed properties of the Lya forest 
in quasar spectra, it follows that this mass must 
exceed about 0.75 keV 10S]. The requirement 
that there be enough small halos to host early 
galaxies to produce the floor in metallicity ob- 
served in the Lya forest systems, and early galax- 
ies and quasars to reionize the universe, probably 
implies a stronger lower limit on the WDM mass 
of at least 1 keV (l09). Simulations pOP do 
show that there will be far fewer small satellite 
halos with AWDM than ACDM. However, as I 
have already mentioned, inclusion of the effects of 
reionization may make the observed numbers of 
satellite galaxies consistent with the predictions 
of ACDM p0fl, in which case AWDM may pre- 



dict too few small satellite galaxies [111]. Lensing 
can be used to look for these subhalos [112,113] 
and may already indicate that th ere a re more of 
them than expected in AWDM fill . Thus it 
appears likely that WDM does not solve all the 
problems it was invoked to solve, and may cre- 
ate new problems. Moreover, even with an initial 



s 



power spectrum truncated on small scales, sim- 
ulations appear to indicate that dark matter ha- 
los nevertheless have density profiles much like 



those in CDM |ll5|j55|,|77j (although doubts have 



been expressed about the reliability of such sim- 
ulations because of numerical relaxation [ 1 16 1 ) . 
But WDM does lead to lower concentration halos 
in better agreement with observed rotation veloc- 
ity curves [ pd7|JTT8| . 

One theoretical direction that does appear 
very much worth investigating is ACDM with a 
tilt n ~ 0.9 in the primordial power spectrum 



Pp(k) oc fc" (ill). Such tACDM cosmology is fa- 
vored by recent measurements of the power spec- 
trum of the Ly a forest |^] and appears to be con- 
sistent with the latest CMB measurements and 
all other available data [119|. Our simple ana- 
lytic model (4^] predicts that the concentration 
of halos in tACDM will be approximately half 
that in LCDM, because the reduced power on 
small scales makes the halos form later. While 
this does not resolve the possible cusp problem, 
it is a step in the right direction which may lessen 
the conflict with galaxy rotation curves. 

8. Conclusion 

The successes of the CDM paradigm are re- 
markable. Except possibly for the density profiles 
at the centers of dwarf and low surface bright- 
ness galaxies, the predictions of ACDM appear 
to be in good agreement with the available ob- 
servations. The disagreements between predic- 
tions and data at galaxy centers appear to occur 
on smaller scales than was once thought. As the 
data improve it is possible that the discrepancies 
on Si kpc scales may ultimately show that CDM 
cannot be the correct theory of structure forma- 
tion. However, ACDM appears to be better than 
any alternative theory that has so far been stud- 
ied, even though these alternative theories have 
additional adjustable parameters. Maybe ACDM 
is even true. 
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